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Abstract 

Radiative double electron capture (RDEC) is a one-step process where two free (or quasi- 
free) target electrons are captured into a bound state of the projectile, e.g. into an empty 
K-shell, and the energy excess is released as a single photon. This process can be treated as 
a time inverse of double photoionization. However, unlike in case of photoionization exper- 
iments, bare ions are used during RDEC observations. Thus, RDEC can be considered as 
the simplest, clean tool for investigation of the electron-electron interaction in the presence 
of electromagnetic fields generated during ion-atom collisions. 

Within this dissertation, the 38 MeV 8+ + C experiment, conducted at Western Michigan 
University using the tandem Van de Graaff accelerator, is discussed and the first experimental 
evidence of the RDEC process is presented. The cross section obtained experimentally is 
compared to the latest theoretical calculations. 
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Abstrakt 

Skorelowany radiacyjny wychwyt dwoch elektronow (RDEC) jest procesem, podczas ktore- 
go dwa swobodne (albo kwaziswobodne) elektrony tarczy wychwytywane sa_ do stanu zwiazane- 
go pocisku (np. nieobsadzonej powloki K), a roznica energii pomiedzy koncowym a poczatko- 
wym stanem eletronow emitowana jest w postaci pojedynczego fotonu. Proces ten mozna 
traktowac jako odwrocenie w czasie podwojnej fotojonizacji. Jednakze, w przeciwieristwie do 
eksperymentow dedykowanych fotojonizacji, do obserwacji RDEC stosuje sie^ jony calkowicie 
pozbawione elektronow, co pozwala na wyeliminowanie tla pochodzaxego od elektronow nie 
biorq,cych bezposrednio udzialu w badanym procesie. RDEC moze bye wi^c traktowany jako 
najprostsze narzedzie do badania oddzialywania elektron-elektron w obecnosci pola elektro- 
magnetycznego generowanego podczas zderzenia. 

Rozprawa ta poswiecona jest procesom atomowym zachodz^cym w zderzeniach 8+ + C 
przy energii 38 MeV podczas eksperymentu przeprowadzonego przy uzyciu alceleratora Van de 
Graaffa w Western Michigan University. Przedstawione zostalo w niej pierwsze doswiadczalne 
potwierdzenie procesu RDEC. Uzyskany eksperymentalnie przekroj czynny zostal porownany 
z wynikami najnowszych przewidywari teoretycznych. 
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Chapter 1 
Introduction 



Since the first observation of the photoelectric effect by Hertz |Her 87] and its explanation 
by Einstein [Ein 05| the interaction between electrons and light has been of considerable 
attention. The fundamental process occurring due to this interaction is photoionization, 
where absorption of a photon of energy fvjj results in the emission of an electron: 

A + hw^A + + e-. (1.1) 

Simple photoionization experiments usually are restricted to neutral atoms, where the influ- 
ence of the electrons, which do not participate in the process directly, cannot be neglected. 
This complicates comparison of the experimental results with theoretical predictions. 

However, based on the principle of detailed balance [Lan 79, Bey 03 1 the photoionization 



can be studied via the time reversed processes, i.e. radiative recombination (RR) and radiative 
electron capture (REC) |Ich 941flch 96|IEic 95aj . During these processes a free (RR) or loosely 
bound (REC) electron is captured to the bound state of the projectile and a photon with 
energy equal to the difference between the final and initial electron states is emitted. Unlike 
single photoionization of multielectron systems, REC has been investigated for bare ion-atom 
interactions |Sto 92| ISto 94 ] and offers clean conditions for exploration of photoionization with 
only one electron, allowing for observation of pure photon-electron interactions. 

During the last thirty years double photoionization has been of considerable interest 
|Dal 94| and references therein]. As a photon typically interacts only with one electron, 
double photoionization is caused by the electron-electron interaction |Smi 89] , However, 
double photoionization studies have been performed mainly for low Z atoms, such as He 
|Ber 931 ITiw 821ICar~8l] . Ne |Sai 921 ISch 93llCar"77j . and Ar |Lab 87llCar"77j . This is due to 
the background contributions from other electrons for high Z systems, which make the subtle 
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electron correlation effects difficult to observe. Fortunately, similar to single photoionization, 
double photoionization can be studied by means of the time inversed process - RDEC, for 
which this background is absent. Radiative double electron capture (RDEC) involves transfer 
of two target electrons into a bound state of the projectile with simultaneous emission of a 
single photon [War 95] IBed 03] . Since bare ions are used during the experiment, RDEC can 
be considered as the simplest, clean tool for investigation of the electron-electron interaction 
|War 95| in the presence of electromagnetic fields generated during ion-atom collision. Thus, 
investigation of the RDEC process can provide crucial information necessary for a proper 
description of electron correlations within atomic systems and provide data required to define 
the wave function of two correlated electrons in the projectile continuum. 

During the last twenty years the RDEC process was addressed not only experimentally 
[War 95|[Bed~03] . but also theoretically |Mir 89| lYaF96| lYak 97] . The calculations were found 
to be in disagreement with the experimental data [Bed 03j and verification of the RDEC pro- 
cess was not possible. The more recent calculations not only explained previous experimental 
results, but also suggested the choice of low energy mid-Z [Z < 35) collision systems for 
observation of RDEC |Mik 04a| IMik 04bl IDru 07] . It is also noted that for these systems 
capture to an excited Is 1 2s 1 state might significantly enhance the process and contribute 
to the observed x-ray spectra [Nef 05] . These calculations provided the main motivation for 
yet another experiment dedicated to the RDEC process. To fully take advantage of the new 
calculations, two collision systems at two different accelerator complexes were chosen: 

• Xe 54+ + C at 20 MeV/u, to be performed at GSI, Darmstadt in Germany, 

• 8+ + C at 2.375 MeV/u, realized by means of the Van de Graaff accelerator at WMU, 
Kalamazoo, MI, USA. 

So far the WMU experiment was carried out. During six months of the experiment prepara- 
tions and data taking, 43 days of beam time were used. At the moment when this thesis is 
being written the GSI beam time is still pending. 

Within this dissertation the 8+ + C at 38 MeV experiment is discussed and the first ex- 
perimental evidence of the RDEC process is presented. The cross section obtained experimen- 
tally is compared with the latest theoretical calculations |Mik 04a. Mi k 04b] INef 05( IDru 07] , 

This thesis begins with an introduction to the most important atomic processes that occur 
during ion- atom collisions. In Chapter special attention is paid to the processes which 
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add to the background for the x-ray spectrum registered during the experiment and formulae 
allowing for estimations of contributions of these processes are suggested. Chapter [3] addresses 
the RDEC process in detail. The history of the experimental approach and the theoretical 
calculations of the RDEC cross section are presented. Additionally this chapter focuses 
on the recent theoretical calculations which were the main motivation for the experiment 
discussed in this dissertation. The goal of the experiment was the observation of x rays emitted 
during collisions of bare oxygen ions with carbon atoms. The x-ray spectra were registered in 
coincidence with ongoing particles which underwent single or double charge exchange. The 
experimental setup which allowed for achieving this goal is presented in Chapter [4] The 
operation principle of a Van de Graaff accelerator is explained and the construction of the 
target chamber, particle spectrometer and x-ray detector are described in detail. Chapter [5] 
is dedicated to data analysis, with a particular focus on processes that may contribute to the 
x-ray spectrum within the RDEC region. Various approaches to estimation of the background 
and calculations of the cross section are discussed. In Chapter [6] the experimentally obtained 
RDEC cross section is compared with the theoretical value and the possible reasons for the 
obtained discrepancy are given. In Chapter [7] results of the Monte Carlo simulations of the 
x-ray spectrum generated during the 8+ + C collisions are compared with the experimental 
results. Finally, in Chapter [8] suggestions for further investigations of the RDEC process are 
given, with indication of necessary improvements of the experimental setup. 
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Chapter 2 

Atomic processes during ion-atom collisions at low 
energy 

Interaction between an incoming ion and a target atom may lead to many different atomic 
processes, such as: 

• ionization, mainly of the target atom, as the electrons are usually less bound to a light 
target than to a partially ionized projectile, 

• electron transfer from the target to the projectile, 

• excitation of both target and projectile states - such states deexcite after the collision 
emitting characteristic x rays. 

Within the following sections the most important processes that were considered compet- 
itive to RDEC for the presented experiment are discussed. 

2.1 Nonradiative electron capture (NREC) 

The Coulomb interaction between the projectile and the target electrons can lead to a 
process called Coulomb capture or nonradiative electron capture (NREC). Here, the energy 
difference between the initial and final state of the electron is converted into the kinetic 
energy of the collision partners. The most convenient and frequently used scaling formula that 
estimates the cross section for nonradiative electron capture is the one given by Schlachter 
|Sch 83] , It is a semiempirical formula which allows for calculation of the (Tnrec as a function 
of the projectile energy for various projectiles with an accuracy better than 30%. 
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The NREC process occurs mainly at the velocity matching condition v ~ v e , where v e is the 
velocity of the captured electron, bound in the target atom. For u > t) e in the nonrelativistic 
approximation, as shown for example in |Eic 07] . the NREC cross section scales as: 

ZfZ 5 

onrec ~ — 19 • (2.1) 
2.2 Radiative electron capture (REC) 

Radiative electron capture (REC) is one of the best known atomic processes observed 
in heavy ion-atom collisions. It was first observed in early seventies of the last century 
|Sch 72 1 IKie 73| ISch 74] and since that time has been intensively studied both experimentally 
|Kan 951 lMok~95l |Spi 79[ ISto 95bl ISto~92l ISto 941 ISto 95al ISto~97bl ISto 97al ISto 981 lTkn~8Tl 



ITan 87] and theoretically |Eic 95al lEic 95bl IHin 871 llch 941 llch 961 ISoh 76j . During this pro- 
cess capture of a single target electron is followed by a photon emission (Fig. 12.1] ). The energy 
Erec of the emitted photon fulfills the energy conservation rule for this process. Thus, it is 
given by: 

Erec = T r + E B - E Bt + (2.2) 

where Eb and Est are the binding energies of the projectile and target, respectively, it 
is the projectile velocity and p the momentum of the electron in the bound state of the 
target. T r = (m e /m p )E is the kinetic energy of the quasifree target electron calculated in the 
projectile's frame of reference. 

The REC line observed during experiments is much broader than the characteristic x- 
ray lines, as can be observed in Fig. 12.21 which is due to the velocity distribution of target 
electrons. This distribution is described by Compton profile Q(p z ) |Big 75 , which gives the 



probability of finding an electron with a given momentum projection p z , where (for ion-atom 
collisions) the z-axis is defined by the projectile velocity. The Compton profile depends on 
the target atomic number Zt and its width increases with increasing Z±. Moreover, the width 
depends on the binding energy and is smaller for loosely bound electrons, than for a tightly 
bound Is electron as shown in Fig. 12.31 

When the binding energy of the target electron is much smaller than T r , the captured 
electron can be treated as quasifree. This means that REC can be described as capture of 
a free electron (radiative recombination - RR), which is the time inverse of photoionization. 
As the cross section for single photoionization can be calculated from the well known formula 
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REC 



Figure 2.1: Radiative electron capture (REC). A target electron is captured into the projectile bound 
state and the energy excess is emitted as a single photon. 



given by Stobbe |Sto 30| . one can use it to calculate the REC cross section via the principle 
of detailed balance. 

Principle of detailed balance describes the relation between the cross sections for direct 
i a i->f) an d time inverse (cr*^) processes [Lan 79, Bey 03 : 

9iPiCi^f(Pi) = 9fPf<?f-n(pf), (2-3) 

where g is the number of possible states given by angular momentum and spin combinations 
and p is the momentum of the particle in the center of mass system describing the size of a 
phase space accessible for the initial (i) and final (/) states. 

Based on Eq. 12.31 and the Stobbe formula for the photoionization cross section, the cross 
section for REC to the projectile K-shell during collision of a bare ion with a hydrogen target 
can be expressed in the form: 



orec = 9-16 



e X p(-4,cot- 1 (l/-)). 10 _ 21 2 



(2.4) 



1 + v 2 ) 1 — exp(— 2iri/) 
where u = Zte 2 /hv is the Sommerfeld parameter of the target K-shell electron and v is the 
projectile velocity. Thus, for fast collisions, the REC cross section scales with energy as: 



Z,Z 5 



greg 



(2.5) 
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125 




180 



E [keV] 

Figure 2.2: Example of the x-ray spectrum registered in coincidence with single electron capture 
during U 92+ + N 2 collisions at 309.7 MeV/u jSwi 00j . 



When this result is compared with Eq. I2.1| one should notice that the radiative electron 
capture dominates for high energy collisions with light targets. 

The angular distribution of the REC photons is given by the angular differential REC 
cross section calculated within the dipole approximation [Se ll 72 1 IKie 73] : 

do rec 3 



, n -o REC sin d. 
ail sir 

Finally, the double differential cross section cP 'a rec I 'dCldE^ can be expressed as: 

d 2 <JREc 1 da rec 



(2.6) 



v dft 



P=P0+Pz 



(2.7) 



where Ss(p z ) is the Compton profile of the target electrons. This formula describes the shape 
of the REC line registered within the x-ray spectrum at a given observation angle. 



2.3 Bremsstrahlung 

When a charged particle penetrates a gaseous or solid target a continuous x-ray spectrum 
is emitted. This spectrum is a result of bremsstrahlung processes occurring in the target 
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Pz 



Figure 2.3: Compton profile of electrons in carbon atom. It can be noticed that the structure of the 
Is line is much broader than that for n = 2 states |Big 75 1 . 




Figure 2.4: Example of the bremsstrahlung process for an electron in the electromagnetic field of an 
ion. 

material, when a charged particle is accelerated (or decelerated) in the Coulomb field of 
the target components. A schematic explanation of this process for an electron in a field 
of an ion is presented in Fig. 12.41 Bremsstrahlung was first observed by Rontgen in 1895 
[Roe 96| IRoe 98] and since that time has been intensively studied |Ish 871 Ush 06| IMir 891 
IChu 811 CTaQgl ILud 98] . 
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RECC 



Figure 2.5: Radiative electron capture to continuum (RECC). Target electron is captured into a 
continuum state of the projectile and a photon is emitted. 

During ion-atom collisions both the ion and ejected electrons may undergo bremsstrahlung 
processes. However, the total power radiated via bremsstrahlung is proportional to 7 4 (when 
dlf/dt _L if) or 7 6 (when dlf/dt\\lf) |Gri 01| . Since E = jmc , where m is the rest mass of 
the moving particle, the total radiated power is proportional to 1/m 4 or 1/m 6 , respectively. 
The above means that electrons lose energy via the bremsstrahlung process much more rapidly 
than heavier charged particles. This is why electron bremsstrahlung dominates over the ion- 
related processes. 

Quasifree electron bremsstrahlung (QFEB), secondary electron bremsstrahlung (SEB), 
atomic bremsstrahlung (AB) and nucleus-nucleus bremsstrahlung (NB) dominate among var- 
ious bremsstrahlung processes that can occur during ion-atom collision. These processes were 
taken into account during data analysis and are more thoroughly discussed in the following 
sections. 



2.3.1 Electron bremsstrahlung 

Radiative electron capture to continuum (RECC), sometimes referred to as quasi-free 
electron bremsstrahlung (QFEB), is a process where the target electron is captured to the 
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projectile continuum, which means it becomes a free electron. Energy conservation in this 
process is fulfilled by a photon emission (Fig. I2.5p . 

The maximum kinetic energy (T r ) of the involved electron, calculated in the projectile 
frame assuming T r S> Est, is given by: 

T r = \m e v 2 = —E, (2.8) 
2 m p 

where v is the velocity of the incoming ion in the laboratory frame (equal to the velocity 
of the captured electron in the projectile reference frame). T r is the maximum energy (in 
the projectile frame of reference) of the photon emitted during the RECC process. As the 
maximum energy of the emitted photons is well defined, the spectrum of the emitted x-rays 
will have an edge at this value. This edge was observed, for example, during collisions of 
carbon ions with Be- and C-targets [Bed 98| . as shown in Fig. 12.61 

Ejected target electrons may scatter in the Coulomb field of other target nuclei, producing 
additional bremsstrahlung. This process is referred to as secondary electron bremsstrahlung 
(SEB). In this case the maximum energy (T m ) of the emitted photons is equal to the maximum 
transfer of the kinetic energy during ion-electron collisions, given by: 

in 

T m = 4T r = A^E. (2.9) 
m p 

Thus, similar to RECC, SEB spectrum has an edge at the photon energy of T m |Ish 87] . 

During the atomic bremsstrahlung process (AB), the projectile excites a target electron to 
a target continuum state. This electron can be recaptured by a target atom with simultaneous 
emission of x rays. The electron can also lose only part of its energy but remain free, in which 
case this process is referred to as radiative ionization (RI). 

It was shown in jlsh 06] that the relative contribution of the above processes strongly varies 
with projectile energy. The theoretical description of bremsstrahlung cross sections given in 
|Ish 06] is in agreement with experimental data. Comparison of experimental data for p + Al 
collisions at 1 and 4 MeV with theoretical calculations are presented in Fig. 12.71 Simple 
scaling formulae describing double differential bremsstrahlung cross sections were proposed in 
jlsh 06] : 

(fay) 2 d 2 a recc _ 7 t ,m e E 
Z2 d(hw)cM ~ tn m p hw h [ ' 

(M 3 d 2 a A B ff flow, r9in 
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Figure 2.6: Spectra observed during the experiment of [Bed 98J for: (a) Be-target, (b) C-target. For 
each energy-target combination two spectra are displayedwith the upper spectrum showing the raw 
data and the lower spectrum with the results after background subtraction. For presentation purposes 
spectra were multiplied by factors; (a) Be-target: 1/20 for 75 MeV/u, 10 for 290 MeV/u; (b) C-target: 
1/8 for 75 MeV/u, 10 for 150 MeV/u, 50 for 290 MeV/u. Dashed line: SEB contribution, solid 
line: RECC (relativistic approximation) + K-REC + SEB. Arrows show the RECC-edge energy T r 
transformed to the laboratory frame. Inset in (a) represents the experimental setup. 



(hiv) 2 d 2 a S EB 



zff( 



m e E , 
m p hio' 



(2.12) 



Z 2 d(huj)dQ 

where fkv denotes photon energy, ao is the Bohr radius and / is a universal function discussed 
extensively in [Ish 06j. The bremsstrahlung processes for protons interacting with various 
targets at a wide range of energies were thoroughly studied for example by Folkmann |Fol 841 
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Figure 2.7: Example of the contribution of various bremsstrahlung processes to the continuous x-ray 
spectrum during pi + Al collisions at 1 and 4 MeV [Ish 06| . It can be noticed that SEB becomes a 
dominating process at higher beam energy. 



IFol 75] . By means of the above formulae, the bremsstrahlung contribution to the experimental 
data can be estimated from the proton data. 



2.3.2 Nucleus-nucleus bremsstrahlung (NB) 

Nucleus-nucleus bremsstrahlung is a consequence of the projectile scattering in the Cou- 
lomb field of the target nuclei. The emitted x-ray spectrum extends up to the projectile 
energy. The differential cross section for NB process can be calculated from the formula given 
by Mokler |Mok 78] : 



d(JNB _ „ Z 2 Z} (Z Zj 
d(huj) ~ (hw)E \A A t 



C = In 



X 



4.3 • 10~ 28 [cm 



A + A t (M 
A t E ' 



(2.13) 
(2.14) 
(2.15) 
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100 



Figure 2.8: Bremsstrahlung processes observed during p + C collisions at 2 MeV. Plot based on Fig. 3 
in |Ish 06] . 

where A, At are projectile and target masses in atomic units, respectively. 

Fig. 12. 81 shows the contribution of RECC, SEB, AB and NB to the x-ray spectra obtained 
during collisions of 2.0 MeV protons with carbon. Upper limits of RECC and SEB x-ray 
spectra can be observed at T r and T m , respectively. The NB cross section is significantly 
lower than those of the electron bremsstrahlung processes. Thus, NB plays a significant role 
only within the x-ray spectrum range above 

2.4 Multielectron capture processes, noncorrelated double ra- 
diative electron capture (DREC) 

During a single ion-atom collision capture of more than one target electron to the projec- 
tile bound state is possible. The simplest example of a noncorrelated capture of two electrons 
is double radiative electron capture (DREC) for which the capture of two electrons is accom- 
panied by the emission of two independent REC photons (Fig. 12. 9p . During this process the 
captured electrons do not interact with each other and the capture of each of them can be 
treated as a separate process. 
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Figure 2.9: Comparison of the DREC and RDEC processes. 



Radiative noncorrelated double capture was theoretically addressed by Meyer hof |M ey 85| . 
In this paper the author calculated the REC cross section as an integral of the process prob- 
ability given as a function of impact parameter. The single electron capture cross section can 
be calculated as an integral of the probability P(b) of an electron capture given as a function 
of the impact parameter b: 

POD 

°rec= \ db2iTbP(b). (2.16) 

J 

In case of REC, P(b) is given by: 



/oo 
dzp(R), 
-oo 



(2.17) 



R being the projectile to target atom distance and p the electron density. As the electron 
density is normalized: 



/•OO /'OO 

/ db2i:b / dzp(R) = Z t , 

JO J -oo 



(2.18) 



the REC cross section for a multielectron target is simply given by Zt(JREc{Zt = !)• 

The same method was applied to noncorrelated double radiative capture. If Po(b) is the 
probability of a single electron capture into the bare ion, and Pi(b) is the corresponding 
probability of electron capture into the H-like ion, the cross section for noncorrelated double 
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1.0 2.0 



Figure 2.10: Single electron loss cross section as given in |Hip 87 . Solid line: PWBA calculations for 
He 2+ impact, dot-dashed line: includes contribution from free electron impact in CBE approximation. 
Symbols: A - Si 8+ , □ - 6 +, V - Si 13+ , o and x - 7+ . 

electron capture can be expressed as |Mey 85 1: 

CTDREC= / db27TbP (b)P 1 (b). (2.19) 
J 

For double electron capture (DREC) one obtains the cross section [M ey 85| : 

(TDREC = 0.13Zl /2 a 2 REC (Z t )a 2 . (2.20) 

The above formula was verified experimentally by Bednarz [Bed 03| , 

When the captured electrons interact with each other during the collision, the process 
is referred to as correlated capture. Radiative double electron capture (RDEC) is the basic 
example of a correlated process and can be treated as a time inverse of double photoionization. 
Thus, due to the electron-electron interaction of the two captured electrons only one photon is 
emitted and its energy is about two times greater than that of a single REC photon (Fig. 12.9]) . 
The RDEC process is discussed in detail in Chapter 
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2.5 Projectile ionization electron loss 

There is a variety of terminology used in the literature when reffering to electron de- 
tachment, which leads to confusion. Here, a nomenclature from [Bom 89] ITan 91| Hip 87 



is 



applied. The term ionization is used when the electron is detached from the target atom, 
while the removal of the electron from the projectile bound state is reffered to as electron loss. 
Consequently, electron loss is the process where an electron is removed from the projectile 
and remains free afterwards: 

A q+ +A t A {q+ '^ + + A t + xe~, (2.21) 

where x is the number of electrons lost by the projectile (q + x < Z). Electron loss processes 
have been extensively studied during the late seventies and eighties for various elements and 
charge states within energy range up to 10 MeV/u |Gra 84| IGra 85| IQls 78| . 

Boman et al. [Bom 89j developed a simple scaling formulae for electron loss cross section. 
The single electron loss cross section for oxygen ions at 1 MeV/u can be estimated as: 

• for q = 5 : 

a\ = (3.27 • l(T 18 )Z t a98 [cm 2 ], (2.22) 

• for q = 6 : 

a\ = (8.83 • l(T 19 )Z t a78 [cm 2 ], (2.23) 

• for q = 7 : 

a\ = (2.22 • 10" 19 )Z°- 33 [cm 2 ], (2.24) 

where q denotes the initial charge state of the ion. It has been also checked by the authors 
that in case of Si 8+ + He collisions at 1.0 MeV/u the ratio of single to double electron loss 
cross sections <t 8 /<7 8 ~ 40. Thus, it can be assumed that the double electron loss process 
can be neglected for the case of the more tightly bound K-shell electrons in 6+ . As can 
be observed in Fig 12.101 the single electron loss cross section does not change significantly 
when the beam energy is increased from 1 to 2 MeV/u. Thus, the scaling formulae given by 
Eqs 12. 22112.241 can be used to estimate the cross sections within this energy range. 



17 



Anna Simon 



Correlated radiative electron capture in ion-atom collisions 



Chapter 3 

Radiative double electron capture (RDEC) 

Radiative double electron capture (RDEC) is a one-step process for which two free (or 
quasifree) target electrons are captured into bound states of the projectile, e.g. into an empty 
K-shell, and the energy excess is released as a single photon (Fig. 12. 9p . This process has to 
be compared with a two-step double radiative electron capture (DREC) during which two 
electrons are captured independently and two single REC photons are emitted. 

While for the DREC process both electrons can be treated separately (see Section 12 .3.2[) , 
in case of the RDEC one has to go beyond the independent electron model. Here, due to 
the electron-electron interaction, transitions of two target electrons into the projectile bound 
states occur with an emission of one photon with an energy about two times greater than 
that of a single DREC photon. 

In general, captured electrons may originate from two different orbitals in the target and 
arrive finally at different final states in the projectile. Thus, the energy of the RDEC photon 
can be expressed as: 

Ebdec « 2T r + + 4 2) - E§ - E% + itfW + (3.1) 

where the indices (1) and (2) correspond to each of the captured electrons. The width of the 
peak is about twice as large as that of the REC line. Roughly, it is determined by the sum of 
Compton profiles of the two active electrons [Mir 89j. 

3.1 Initial experiments dedicated to RDEC 

The first experiment dedicated to RDEC was performed at GSI in 1994 with 11.4 MeV/u 
Ar 18+ ions from the UNILAC impinging upon a carbon foil. A detailed description of this 
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Figure 3.1: Typical x-ray spectrum obtained during argon experiment [W ar 95) . 



experiment is given in [War 95] . A typical spectrum obtained during that experiment is 
presented in Fig. 13.11 As shown in this figure, no significant line structure related to the 
RDEC process was observed. However, the number of counts collected in the expected RDEC 
energy window provided an upper limit of the RDEC cross section of about 5.2 mb. A rough 
theoretical estimate of the cross section ratio (Trdec I '&REC was a l so suggested, based on the 
principle of detailed balance and considering REC as a time reversal of photoionization. The 
REC and RDEC cross sections can be written as: 

Huj 



prec{(3) = Zt 



api(hw) 



VRDEC(P) = FZ t (Z t - 1) 



a DPI (hw') 



(3.2) 



(3.3) 



2'yf3mc 2 1 

where a pi and a dpi are the cross sections for single and double photoionization, respectively. 
The factor F {F < 1) describes the phase space fraction of double photoionization accessible 
for the RDEC process. Thus, the ctrdec/^REC ratio can be expressed in terms of single and 
double photoionization cross sections as [War 95j: 



R 



crdec 



F(Z, 



f u/_\ 2 f a DPI {hijj') 



prec 

or, as in case of RDEC tkj' ~ 2Huj: 

CRD EC 



R 



&REC 



F(Z t - 1) 



\2uj J \ api(Huj) 

apj(2hw) aDPi(2fwj) 
api(huj) api(2huj) 



(3.4) 



(3.5) 
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Figure 3.2: Typical experimental x-ray spectrum obtained for uranium ions [Bed 03] . The Gaussian 
solid line shows the expected RDEC peak, which should be observed according to Yakhontov et al. 
| Yak 96llYak~97] . 

where the values of the expresssions: a pi(2hui) / 'a pi{ftw) and aopi{2hw) / 'apj(2huj) can be 
easily estimated from [Fan 591 1 Amu 75] , The experiment |War 95] provided an upper limit 
for R of 3.1 • 1(T 6 . 

This experiment stimulated theoretical treatment of the RDEC process |Yak 96 1 lYak 97] . 
In these papers the authors presented nonrelativistic calculations of the RDEC process adapted 
to the kinematics and the energy range of the Ar 18+ + C experiment. The calculations gave, 
for this particular collision system (Ar 18+ + C at 11.4 MeV/u), the RDEC to REC cross 
section ratio R of 3.6 • 10 -6 , which is close to the experimental upper limit. 

Moreover, these calculations predicted a strong enhancement of RDEC during heavy ion- 
atom collisions at relativistic energies |Yak 97] . These calculations were tested during the 
second experiment dedicated to RDEC. Here, bare uranium ions at an energy of 297 MeV/u 
collided with an Ar target at the ESR storage ring of the GSI facility [B ed 03] . This experiment 
showed that for the collision system under consideration the RDEC cross section is certainly 
at least three orders of magnitude smaller than the theoretical prediction |Yak 96 1 lYak 97] , 
Fig. l3.2l shows a spectrum obtained during the experiment. Again, no significant line structure 
which could be assigned to the RDEC process was observed. The Gaussian line in Fig. 13.21 
shown within the RDEC region of the spectrum represents the shape of the RDEC line which 
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should be observed in the spectrum, if the theoretical calculations [Yak 96, lYak 97| were 
reliable. This experiment also provided only an upper limit for the RDEC cross section value 
of about 10 mb. 



3.2 Recent theoretical approach to RDEC 

In order to explain the disagreement between the uranium experiment [Bed 03j and the 
theoretical treatment of RDEC |Yak 96|lYak~97j . a new theoretical approach for the correlated 
double electron capture into the K-shell of bare ions was proposed |Mik 04a"l IMik 04b"l INef 05] . 
Here, a brief description of this RDEC treatment is given with the notation used in the original 
papers. Indices (1) and (2) correspond to REC and RDEC, respectively, and natural units 
(h = c = 1) are used throughout the text. 

All the electrons involved in the process were considered as nonrelativistic and the energy 
oj of the emitted photon was limited by 1<ik < oj m, where I^k is the threshold energy 
for double photoionization of the projectile K-shell and m is the electron mass. In such case 
the Coulomb parameter (aZ, a denotes the fine structure constant) is small {aZ <C 1) and 
perturbation theory with respect to the electron-electron interaction can be used. 

In the reference frame of the incident ion the probability dW for double electron capture 
into the K-shell of bare ion with the emission of a single photon per unit time is given by 
|Mik 04a| : 

9tt die 

dW = —\Af^5(2Ep -u- I 2K ), (3.6) 

where Ep is the one-electron energy within the initial continuum state, u) = | k \ = k is the 
energy of the emitted photon and I2K = 21 with I = rj 2 /2m, the Coulomb potential for single 
ionization and rj = maZ the characteristic momentum of the K-shell electron, and V is a 
normalization factor. Summation over all polarizations of the emitted photon is assumed in 
Eq. 13.61 and the delta function ensures energy conservation. The amplitude A was obtained 
from that for the double K-shell photoionization. A detailed description of this approach is 
given in [Mik 04aJ. 

Dividing Eq. 13.61 by the current flux of the incident target electrons j = v/V, where 
v = p/m is the absolute value of the initial velocity of the incident electrons before collision 
with the ion, one obtains the effective differential cross section: 

d*V=2n^\A\^ (3.7) 
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Figure 3.3: Universal function Q calculated as a function of the dhnensionless variable £ [Mi k 04a] , 



which defines the angular distribution of the RDEC photons emitted into an element of a 
solid angle dSl^. 

For the energy regime assumed in these calculations, it was possible to calculate the total 
cross sections within the electric dipole approximation. For collisions of heavy ions with light 
target atoms the total cross section for radiative double electron capture (RDEC) into the 
K-shell of the ion is given by: 



oi9 73 
r (2) = 2 z t 

3vrZ 5 



Q(0, 



(3.8) 



where ( = jj/p is a dimensionless parameter, a$ = o 3 Oq and ao denotes the Bohr radius. 



Q is a universal function of £, which can be obtained by numerical integration (Fig. I3.3|) . 
£ ~ 1 corresponds to the near-threshold domain, where the K-shell photoeffect reaches its 
maximum. For slow collisions (£ S> 1) the RDEC cross section increases, while in case of fast 
collisions it decreases rapidly. Moreover, it has to be pointed out that the RDEC cross section 
rapidly drops with the projectile atomic number (~ Z -5 ) and increases significantly for low 
energy collisions. 

Another value calculated in (Mik 04a| is the cross section ratio (R = cr^/crW). The 
REC cross section can be expressed in terms of the photoionization cross section (<7pj) via 
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Figure 3.4: Universal quantity Q/H calculated as a function of the dimensionless variable £ |Mik 04a] . 



the principle of detailed balance. As api is known analytically (Stobbe formula |Sto 30| ). for 
radiative electron capture into the K-shell of the projectile one obtains: 



olO 

cj(1) = -y^oZtHit), 



H(0 



£ 2 exp(-4£ cot -1 ^ 
if l-exp(-27rf) ^ 



(3.9) 



(3.10) 



where e 7 =o;// = 2(l + £ 2 )is the dimensionless photon energy. Then the ratio R is given 
by: 



R 



(3-11) 



TT 3 Z 5 H{^' 

The function Q(£)/H(£) is presented in Fig. 13.41 

These calculations are in disagreement with the previous relativistic approach [Yak 97] . 
which was not able to explain the existing experimental data [Bed 03] . As shown in [Mi k 04a] . 
the enhancement of the wave function for the relativistic systems was calculated incorrectly 
by Yakhontov [Yak 97] and even the corrected value, which is 3 orders of magnitude smaller 
|Mik 04a] . was not confirmed by the experiment [Bed 03] , Therefore, the enhancement of the 
RDEC cross section for relativistic systems |Yak 97] seems to be absent. This is similar to 
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Table 3.1: Comparison of experimentally obtained RDEC cross sections [War 9~5l IBed 03| and the 
calculated values given in [Yak 97] and [Mik 04a] , 



z 


E [MeV/u] 




z t 




a( 2 ) [mb] 




Ref. [Mik 04a| 


Ref. |Yak97| 


experiment 


18 


11.4 


0.840 


6 


3.2 


1.85 


<5.2 [War 95| 


92 


297 


0.841 


18 


2.540~ 2 


5000 


<10 [Bed 03) 





1 ■ ■ ■■■■■■ 

1 10 

f 

f2l (2) 

Figure 3.5: The ratio of the RDEC cross sections to the excited (cr^ig) an d ground (cr lls ) projectile 
states as a function of adiabacity parameter £ [Nef 05] . 

the behavior of the cross section for the REC process, where the cross section decreases when 
the projectile energy increases. 

However, it has to be emphasized that the current estimate [Mik 04a] of gives values 
closer to the experimentally obtained upper limits for both the nonrelativistic case (Ar 18 + + C, 
[War 95] ) and the relativistic one (U 92+ + Ar, [Bed 03j) (see Table EE), which suggests that 
[Mik 04a"] is so far the most reliable theoretical description of RDEC. 

In contradiction to predictions given in [Yak 97j . the new calculations show that the RDEC 
cross section strongly depends on the target atomic number and electron density. One can 
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Table 3.2: The REC (<rW), RDEC {a {2 ^) and DREC (a^ 2 ^) cross sections and their ratios as given 
in [Dru 07j . 



z 




E [MeV/u] 


z t 


<7« [kb] 


[mb] 


a (2 ^ [mb] 


CT (2,7) /(J (1) 


(T (2,7)/ (T (2,2 7 ) 




0.84 


11.4 




0.36 


1.5 


3.2 


8.9-10- 6 


2.1 


18 


0.20 


646 


6 


1.5-10~ 3 


2.6-10~ 5 


LO-IO" 6 


6.7-10~ 10 


3.2-10~ 2 




0.10 


804 




6.4-10~ 5 


4.7-10~ 8 


1.6-10- 10 


4.0-10~ 12 


3.4-10- 3 




0.84 


5.1 




0.36 


1.5 


24 


6.7-10- 5 


16 


12 


0.20 


287 


6 


1.5-10~ 3 


2.6-10~ 5 


7.6-10- 6 


5.1-10~ 9 


0.29 




0.10 


357 




6.4-10" 5 


4.7-10" 8 


1.2-10" 9 


1.9-1Q- 11 


2.6-KT 2 



expect much larger values of in case of slow collisions of multicharged ions with a solid 
state target with low atomic number Zt [Mik 04a"] . As the orbital velocity of the target valence 
electrons is much smaller than that of the projectile, they can be considered as quasifree in the 
projectile's frame of reference. In this reference frame these electrons appear as an electron 
beam with velocity v and concentration n e = KptN^/Mt, where k is the number of valence 
electrons, Na is Avogadro's number and pt and Mj are the density and molar mass of the 
target, respectively. Hence, by substituting V = n~ l in Eq. 13.71 q( 2 )can be expressed as: 

° {2) = (n e a s )^Q(0. (3.12) 

In addition, the correlated double electron capture into the ls2s state increases the cross 
section for the RDEC process |Nef 05] , As shown in Fig. 13. 5| the ratio of the cross section 
for RDEC to the ls2s state, <JgL S , to the cross section for RDEC to the Is 2 ground state, 
Cjio! is strongly dependent on the £ value. As can be seen from Fig. 13.51 for £ S> 1 (i.e. slow 
collisions) the cross section for electron capture to the ls2s state can greatly exceed that for 
Is 2 state capture. 

Recently, the calculations of Nefiodov and Mikhailov were continued by Drukarev |Dru 07] , 
who again addressed the high energy nonrelativistic limit (£ ^> 1) of the RDEC process. As 
previously found, a strong energy dependence of the cross section was shown and the RDEC 
probability was compared with the one for noncorrelated capture. Obtained values of the 
REC, RDEC and DREC cross sections for Ar 18+ + C and Mg 12+ + C for various projectile 
energies are given in Table I3~2"1 

This theory [Mi k 04a| IMik 04bl INef 05 \ IDru 07] suggests that the best systems for obser- 
vation of the RDEC process are low energy collisions of mid-Z ions with light solid targets. 
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This theory was a motivation for the next RDEC experiment and a reason for the choice of 
the conditions of the experiment presented in this thesis. 
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Experimental setup at Western Michigan University 
4.1 Van de Graaff accelerator 

The Van de Graaff accelerator is an electrostatic generator which uses a moving belt to 
accumulate very high, electrostatically stable voltage on a hollow metal sphere [Gra 47J. This 
type of generator was developed by Robert J. Van de Graaff at Princeton University. The first 
model was demonstrated in October 1929 and in 1931 a version able to produce a potential 
difference of 1 MV was described |Gra 31] . 

A simple Van de Graaff generator is presented in Fig. 14.11 A belt of dielectric material 
runs over two rollers, one of which is surrounded by a hollow metal sphere the high voltage 
terminal. Two electrodes, an upper and a lower one, are placed next to each roller. The upper 
electrode is connected to the sphere, while a high DC potential (with respect to the ground 
potential) is applied to the lower one, a positive potential in the example. 

Due to the strong electric field the air around the lower electrode is ionized and the positive 
ions are repelled from the electrode and accumulated on the belt. Then they are transported 
towards the upper electrode which collects the charges from the belt and transports them 
onto the spherical collecting electrode. The potential of the HV electrode increases until the 
speed of its charging equals the speed of discharging. The maximum potential obtained on 
the HV electrode depends on the radius of the sphere and insulating properties of the gases 
surrounding it. SFg or a mixture of N2 and CO2 under a pressure even up to 20 bar are 
usually used |Hin 97] . The value of the terminal voltage in Van de Graaff accelerators may 
reach up to 15-20 MV |Edw 931 |Bey 03| . 

If a source of positive ions is placed close to the high voltage terminal, as in Fig. 14. 1| the 
ions are repelled by the positive charge of the terminal electrode and thus accelerated towards 
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Figure 4.1: Schematic view of a classical Van de Graaff accelerator: (1) lower roller, (2) upper roller. 
(3) charging electrode, (4) electrode collecting positive charge, (5) voltage generator, (6) spherical 
electrode (high voltage terminal), (7) ion source, (8) extracted ion beam. 

the ground potential. The final kinetic energy of the ions depends on their charge state q and 
is proportional to the terminal voltage Vt er minaV- 

E = qVterminal- (4.1) 

In modern ion accelerators with a Van de Graaff generator, electrodes located at entry and 
exit of the vacuum tube are grounded and the high-voltage terminal is located at the middle 
of the tube, as shown in Fig. 14.21 [Hin 97| [Wed 99J. A source of negative ions is placed at the 
entrance of the tube and produced ions, usually singly charged, are accelerated within the 
tube towards the high- voltage terminal, where two or more electrons are removed from each 
ion as it passes through a stripping foil. The charge state of the ion changes from negative 
to positive and the ion is repelled from the terminal and accelerated towards the grounded 
exit of the tube. Compared to Van de Graaff accelerators of the ordinary type, by means of 
tandem Van de Graaff accelerators higher particle energies can be obtained since the potential 
difference is used for the acceleration twice. Thus, in this case the final kinetic energy can be 
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Figure 4.2: Schematic view of a tandem Van de Graaff accelerator: A negative ion entering the 
accelerator (A - ) is accelerated by the high terminal voltage. Some of its electrons are removed while 
the ion passes through the stripping foil. The positive ion (A +9 ) is repelled by the high voltage 
terminal, thus causing additional acceleration. 

calculated as: 

E= (q + l)Vterminal, (4.2) 

where q is the ion charge state after passing through the stripping foil. 

The experiment described in this dissertation was performed at Western Michigan Univer- 
sity using 6 MV tandem Van de Graaff accelerator. The WMU accelerator was built by the 
High Voltage Engineering Corporation, the company founded by Robert Van de Graaff. The 
construction of the accelerator allows for obtaining stable beams of bare ions of all elements 
up to 9 F with a total kinetic energy up to 40 MeV. 



4.2 Beam line setup at Western Michigan University 

A schematic view of the WMU accelerator beam line is presented in Fig. 14.31 The acceler- 
ated beam passes through a 90° analyzing magnet which allows for choosing the appropriate 
ion charge state. At this point, the final energy of the beam is defined accordingly to Eq. 14.21 
Then the beam passes through a post stripper followed by a switching magnet which directs 
the desired charge state towards the experimental area. For the presented experiment a beam 
of 6+ was extracted from the accelerator operating at the terminal voltage of 5.43 MV, 
which produced a beam of energy equal to 38 MeV. Then the beam traversed through a 20 
jiig/cm 2 carbon stripper foil, where bare and H-like oxygen ions were produced. Simply, by 
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Figure 4.3: Schematic view of the WMU van de Graaff accelerator facility |Kay|. 



changing the magnetic field of the switching magnet one could choose the necessary charge 
state. When a proton beam was accelerated, the post-stripper was removed from the beam 
line. 

The experimental beam line farthest left, when looking along beam direction, was used 
during this experiment. There, an exclusively designed chamber for a solid target was placed, 
which not only allowed for mounting up to four films but also provided a simple mechanism 
for target rotation. This was necessary for optimization of the target position during the 
experiment. During data acquisition the target film was positioned at 45° to the beam di- 
rection, facing the x-ray detector as shown in Fig. 14.41 This setup ensured a direct detection 
of emitted photons, as they did not traverse through the foil, so the unnecessary energy loss 
was avoided. It also allowed for usage of the whole active area of the x-ray detector, which 
was not covered by the aluminum frame of the target holder. The target foils used during 
the experiment were a few //g/cm 2 thick which corresponds to an areal density on the order 
of 10 17 particles/cm 2 . 
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Figure 4.4: The experimental target chamber in 1:1 scale. 



The target chamber was designed in a way that minimizes the distance between detector 
window and target center. The total crystal-target distance achieved was about 25 mm, which 
gives a detection solid angle of AQ = 0.044(1) sr. 

Emitted x-rays were registered by an ORTEC single crystal Si(Li) detector placed perpen- 
dicular to the beam direction. The crystal of 6 mm diameter and 3 mm thickness, together 
with a 7.5 |im Be- window, gave a detection efficiency in the energy range 2-4 keV better than 
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Figure 4.6: Experimental setup. 



/o (Fig. 14. 5p . The detector was energy calibrated with a standard 55 Fe calibration source. 
Calibration was frequently repeated throughout the experiment in order to check the stability 
of the data acquisition system. 

Along the beam direction, a set of two collimators was placed in front of the target 
chamber. The distance between collimators was about 2 m. Collimator apertures of 2 and 
3 mm were to ensure a good beam collimation. An additional collimator between the target 
and magnet prevented scattered ions from entering the spectrometer and generating false 
coincidences (see Fig. 14. 6|) . 



32 



Anna Simon 



Correlated radiative electron capture in ion-atom collisions 



The target chamber was followed by a magnetic spectrometer. The magnetic field of 
the dipole magnet separated the final charge states of the ions and directed them towards 
four surface barrier detectors. The primary beam was registered by a Faraday cup. Surface 
barrier detectors counted ions with charge states equal to q — 2, q — 1, q + 1 and q + 2, where 
q is the charge state of the primary beam. Both the magnetic field of the spectrometer and 
the position of the surface barrier detectors were adjustable and created a versatile system, 
which could be used for various beam charge states and energies. A schematic drawing of the 
experimental setup is presented in Fig. 14.61 

4.3 Data acquisition system 

The data acquisition system was designed for registration of x-rays as well as particles 
with final charge states of q — 2, q — 1, q + 1 and q + 2 in a way that allowed for a software 
analysis of photon-particle coincidences. A schematic drawing of the electronics is presented 
in Fig. 14.81 For the clearness of presentation only a setup for one particle detector is drawn. 

Signals from the particle detectors were first amplified by ORTEC 474 timing filter am- 
plifiers (TFA), converted to a logic signal by constant fraction discriminators (ORTEC 463 
or 473) and finally used as the STOP signals for time to amplitude converters (TAC, OR- 
TEC 566). 

A signal from the x-ray detector preamplifier was amplified and split (into two unipolar 
and one bipolar signals) by the Tennelec 244 spectroscopy amplifier. Additionally, a signal 
from the backside input of the amplifier was used (it is an unchanged signal from the front 
input) as an input for the TFA (ORTEC 474). The unipolar signals were processed by linear 
gate stretchers (LGS, ORTEC 542), while the bipolar signal was analyzed by a timing single 
channel analyzer and then converted to a logic signal by the delay generator (Phillips Sci 794), 
so that it could be used as a strobe signal for TAC, LGS modules and the ADC. The same 
signal was also used as a master trigger for the computer. Additionally, the TFA signal was 
used as a gate for one of the LGSs. This resulted in registration of only x rays which gave the 
START signal for the TAC. The 'non-gated' LGS served as a check of the x rays lost by the 
TFA. It was seen from the data that about 25% of the x-rays were not accepted by the TFA. 

In Fig. 14. 81 a signal from the 'gated' LGS is called 'fast', while the one from the 'non-gated' 
LGS is referred to as 'slow'. The x-ray timing signal was used as a START signal for TACs. 
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Figure 4.7: Example of a time spectrum registered during the experiment. The arrow indicates the 
width of a time window for true coincidences (calibration 2 ns/channel). 



All the TAC and LGS outputs were registered by the ORTEC 811 ADC and the data 
were written into SpecTCL . evt list mode files and converted afterward into . ROOT files for 
analysis purposes. The CERN-ROOT software was used for data analysis. 

The construction of the data acquisition system allowed for registration of the time differ- 
ence between photon and particle detection with a resolution of 86(4) ns (FWHM of the time 
peak, Fig. 14. 7p . This information was used for determination of the true coincidences between 
registered photons and ions. An example of a time spectrum registered during the experiment 
is presented in Fig. 14.71 The peak associated with true coincidences is clearly visible and the 
time window set for data analysis is marked with an arrow. The remaining parts of the time 
spectrum include random coincidences that were, after normalisation, subtracted from the 
coincidence spectrum. 
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Figure 4.8: Scheme of the data acquisition system. 
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Chapter 5 
Data analysis 



The main goal of the experiment was observation of the x rays generated by bare oxygen 
ions (0 8+ ) impinging on carbon at an energy of 38 MeV. This gives a projectile velocity 
of 9.7 a.u., which corresponds to the adiabacity parameter £o of 0.82. As can be seen in 
Fig. 15.11 (a), the single x-ray spectrum registered during 8+ + C collisions is dominated by 
the Al K-a line, which is produced by scattered ions hitting the aluminum target holder. In 
order to establish the shape of this line a run without the carbon foil was performed. The 
resulting spectrum, which contains only the Al K-a line, is represented by the dashed line in 
Fig. 15.11 (a). In Fig. 15.11 (b) the 8+ + C spectrum after subtraction of the Al K-a line is 
shown. 

A run with 38 MeV 7+ ions was made to check if the structure of the x-ray spectra 
in the RDEC region changes when one electron in the projectile K-shell is present. As the 
experiment was aimed at observation of RDEC mainly to the ground state (Is 2 ), presence of 
a Is electron in the H-like ion should block this process. Thus, a change in the RDEC range 
of the x-ray spectrum should be registered. As shown in Fig. 15.11 (c), the structure in the 
RDEC region is different for the 7+ ions. 

As an additional test of the experimental conditions, PIXE analysis of the target foil was 
performed with 2.375 MeV protons, i.e. for the same collision velocity as for 8+ ions. During 
this process the impinging proton beam excites atomic states of the atoms within the material, 
which then deexcite emitting characteristic x rays. Intensities of the lines allow for estimation 
of the amount of the impurities in the material. Here, PIXE analysis was performed to check 
if any impurities that might produce x rays in the RDEC range are present in the foil. As 
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can be seen in Fig. 15.11 (d), no structure in the RDEC region was observed, nor was there 
evidence for REC around the photon energy of 2 keV. 

Coincidence spectra obtained for bare 8+ ions capturing one or two electrons are pre- 
sented in Fig. 15.21 These spectra were obtained using the coincidence condition on the ap- 
propriate time spectra as described in Section 14.31 The contribution of random events was 
subtracted. In both spectra, a well separated structure in the RDEC region is evident. The 
REC counts associated with ions which captured two electrons come from multiple capture 
processes in the target with at least one radiative capture (REC). This could not be avoided 
due to a very high cross section for nonradiative electron capture, which is of the order of 
0.2 Mb, as estimated according to the commonly used scaling formula |Sch 83] , Another origin 
of these REC counts is the DREC process, which, due to a cross section of 48 b (Eq. I2.20p . is 
a source of about 10% of the total number of counts in the REC energy range. 

The following subsections of this chapter will address the data analysis in detail. 

5.1 PIXE analysis of the target material 

To check if there is any statistically significant structure in the RDEC range of x-ray 
spectrum obtained for the proton beam, the method described in Appendix lAl was used. The 
value of a was set to 0.05, which means that the probability that any structure in the tested 
range of the spectrum is due to statistical fluctuations (Hq hypothesis) is equal to 0.95. A 
very crucial step of this analysis was a proper estimation of the background in the vicinity 
of the RDEC range. As can be seen in Fig. 15.31 the background shape is not very smooth 
and a standard linear fit of the background, implementing the least squares method, was very 
sensitive to the choice of the fitting region. That is why the following method of estimation 
of the background parameters was used. A region on both sides of the RDEC range (marked 
with dashed lines in Fig. 15. 3[) . without the Al K-a line, was integrated and the result was 
divided by the number of channels in the chosen range of the spectrum. This gave an average 
number, N c h, of the background counts per channel, N c h = 6.40(23). This number multiplied 
by the number of channels in the RDEC region of the x-ray spectrum gave the total number 
of the background counts of b = 895(32). The total number, N, of counts registered within 
the RDEC region of the spectrum is equal to 903(30). 
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Table 5.1: Results of a X 2 test of the RDEC range of the proton induced spectrum. 



N ch 


N 


b 


T 


1 — cut 


6.40(23) 


903(30) 


895(32) 


0.372 


0.54 



The above values of N and b determined the value of the statistical variable T (Eq. IA.4|) 
to be equal 0.372. This gives 1 — cxt = 0.54, which means that the Hq hypothesis cannot be 
rejected. In other words, there is no evidence for any peak structure, that could be a result 
of any physical process, in the RDEC range of the proton induced spectrum. Thus, it can 
be assumed that there are no target impurities that could emit photons in the RDEC energy 
range. Values of all the relevant parameters are given in Table 15.11 

5.2 Projectile K- and L-shell electron loss 

It could be already noticed that the events in the RDEC range are more frequent in 
the spectrum associated with the single charge exchange rather than with the double charge 
exchange channel (Fig. 15 . 2 1) . Similarly, most of the REC photons are observed in the singles 
spectrum and not in the single charge exchange channel (compare Figs 15.11 (b) and 15.21 (b) ) . It 
is likely that after the observed capture process the ion may lose one of the captured electrons. 
This is due to the fact that, even in a very thin carbon foil, the electron loss probability is 
very high for the weakly bound ionic systems that are formed during the collision. Electron 
loss cross sections obtained by interpolation of relevant data (see Chapter [2]) are given in 
Table I5T21 These cross sections are in fair agreement with the data of Shima et al. [Shi 92] . 
where the bare and H-like ions are indicated as the most populated charge states observed 
when a 38 MeV oxygen beam traverses through a carbon foil. Together, they account for 
more than 80% of the final charge states, with 50% in the H-like (7+) state. For the system 
investigated, the electron loss cross sections are at least 3 orders of magnitude greater than 
the REC cross section. This means that there is a very high probability for an ion, after the 
capture process, to undergo another collision, which is most likely to be accompanied by the 
electron loss process. 

For the system investigated, the cross section for the removal of the L-shell electron is 
about an order of magnitude greater than that for the K-shell electron (see Table [572]) . Thus, 
in case of double capture to the ls 1 2s 1 state, the 2s electron is promptly removed, while in 



38 



Anna Simon 



Correlated radiative electron capture in ion-atom collisions 



Table 5.2: Electron loss cross sections for oxygen ions at 38 MeV estimated from the data presented 
in |Bom 89llTan~9Tl|Hip 87| . 



Process Cross section [Mb] 



L-shell ionization 






Q5+ _ 


-)■ 6 + 


19.0 








Q6+ _ 


■> 7 + 


3.6 


K-shell ionization 


< 




7 + - 


-> 8 + 


0.4 








Q6+ _ 


-> 8+ 


0.1 



case of capture to Is 2 the final charge state of the ion is more likely to remain unchanged. 
Therefore, one would expect most of the photons originating from RDEC to the projectile 
excited state to show up in the single charge exchange channel, while the capture to the 
ground state will be less affected by the K-shell electron loss process. This can be observed in 
Fig. 15.21 where the ls 1 2s 1 peak is clearly visible in the q — 1 coincidence spectrum, while it is 
almost absent in the q — 2 channel, compared with the structure associated with the capture 
to the Is 2 state which is still visible in the q — 2 channel. 

The cross section for double electron loss (0 6+ — > 8+ ), even though it is much smaller 
than that for single electron loss, was also considered. Due to this process, in some cases, both 
captured electrons could be removed from the ion and the RDEC events might be observed 
only in the singles spectrum. However, after subtraction of all the background contributions 
(Al K-a, REC, bremsstrahlung) the remaining number of counts within the RDEC range of 
the singles spectra was consistent with the sum of all the counts from the q—1 and q — 1 
coincidence channels. This means that these two channels include all the registered events 
that could be associated with the RDEC process and the sum of the RDEC counts from both 
channels should be included in the cross section estimation. In case of REC, in order not to 
lose any of the registered events, the total number of REC photons was obtained from the 
singles spectra. 

5.3 Background processes 

For such a subtle effect as the RDEC process, estimation of the background shape in the 
RDEC range of the x-ray spectrum is of a great importance. The most significant processes 
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occurring during collisions of bare ions with solid target that can contribute to the observed 
x-ray spectrum, are: 

• radiative electron capture to continuum (RECC), 

• secondary electron bremsstrahlung (SEB), 

• atomic bremsstrahlung (AB), 

• nucleus-nucleus bremsstrahlung (NB), 

• the high energy tail of the single REC profile caused by the momentum distribution of 
the target electrons (Compton profile). 

Under the experimental conditions, the high energy limits for the RECC (T r ) and the 
SEB (T m ) processes are 1.3 keV and 5.2 keV, respectively (see Eqs 12.81 and 12 . 9[) . Therefore, 
RECC will not contribute to the RDEC background, as it appears in the low photon energy 
range. In addition, RECC will be suppressed by the absorption in the Be-window. The SEB 
process might be a significant background contribution and will be discussed further. 

There are three bremsstrahlung contributions in the x-ray spectra within the RDEC 
photon energy range: SEB, AB and NB. Secondary electron bremsstrahlung and atomic 
bremsstrahlung can be estimated by simple scaling from the data obtained by Folkmann for 
protons of similar velocity to the oxygen ions used during the experiment. The x-ray spectrum 
for p + C collisions given in |Ish 06] can be easily transformed by means of appropriate scaling 
formulae (Eqs 12.10112.12]) . As can be seen in Fig. 12.81 this x-ray spectrum is dominated by 
electron bremsstrahlung (SEB + AB) for photon energies up to 10 keV, which is far beyond 
the region of interest (2.8-4.2 keV). According to the formula given by Eqs l2.13ti2.15l the cross 
section for the NB process is of the order of 1 mb/keV sr in the RDEC range, which is a few 
orders of magnitude smaller than the cross sections for any of the electron bremsstrahlung 
processes. Thus, the contribution from NB can be neglected in the discussed case. 

Another contribution that should be taken into account is the high energy tail of REC. 
The observed REC structure originates not only from single radiative electron capture, but 
can also be associated with the noncorrelated double electron capture (DREC). However, 
estimation based on Eq. 12.201 shows that the DREC cross section is an order of magnitude 
smaller than that for a single REC process, as shown in Table [5T31 
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Table 5.3: Total cross sections for the background processes that were taken into account during data 
analysis. 



Process 


Cross section [kb] 


Radiative electron capture (REC) 


0.512 


Nonradiative electron capture (NREC) 


240.0 


Double radiative electron capture (DREC) 


0.048 



All these background contributions are shown in Fig. 15.41 In the expected RDEC energy 
region (2.8-4.2 keV) the broad REC spectrum significantly exceeds the bremsstrahlung con- 
tribution (see inset in Fig. 15. 4p . As a consequence, the RDEC line will be placed mainly on 
the REC tail. 

The cross section for background contributions was also corrected for the detection effi- 
ciency of the Si(Li) detector. Both crystal efficiency and Be-window absorption were taken 
into account by using the data shown in Fig. 14.51 The shape of the background, after effi- 
ciency correction, is also shown in Fig. 15.41 It can be noticed that the spectrum is completely 
dominated by the tail of the single REC line and the contribution of the bremsstrahlung is 
almost negligible. Additionally, the low energy photons (energy not exceeding 0.5 keV) are 
completely absorbed in the beryllium window of the x-ray detector and the curve, representing 
the sum of all background contributions corrected for the detector efficiency, drops almost to 
zero at low photon energies. 

5.4 Pile-up of single REC photons and its contribution to the 
RDEC energy range of the spectrum 

As the energy of a RDEC photon is about twice as large as the energy of a single REC 
photon, it might be impossible to distinguish a real RDEC event from a situation when 
two REC photons are simultaneously registered by the x-ray detector (pile- up effect). This 
process may produce an additional background within the RDEC range of the x-ray spectrum. 
Simultaneous detection of two REC photons may occur in three situations: 

• both photons emitted during the DREC process propagate in the same direction and 
are both registered by the x-ray detector, 
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Table 5.4: Probabilities and count rates of the processes that might contribute to the x-ray spectrum 
in the RDEC range. For more information see text. 

Process Probability Count rate [events/s] 

RDEC 5.3-10" 11 

DREC 6.04CT 11 

double REC collisions 5.2-10 -14 
sequential REC 



1.6-10" 6 
1.8-HT 6 

i.6-i(r 9 
2.2-icr lc 



• the ion undergoes multiple collisions during which at least twice the REC process occurrs 
and both photons are emitted towards the x-ray detector, 

• if the beam intensity is high enough, there is a chance that two REC photons emitted 
by two sequential ions are simultaneously detected by the x-ray detector. 

If the cross section for RDEC is of the order of 0.1 b |Nef 09| and the given geometry of 
the experiment is taken into account, one obtains the probability of observation of a RDEC 
photon during an ion-atom collision equal to 5.3 • 10~ n . If the beam intensity is equal to 
Imedium = 3 • 10 4 ions/s, which is the mean value obtained during the experiment, the RDEC 
count rate will be of about 1.6 • 1CP 6 events/s for the geometry applied (Table l5T4|) . 

For a low-Z system, as the one used during the experiment, the cross section for the 
noncorrelated double electron capture (DREC) is rather high (see Table 153]) when compared 
to the RDEC cross section. Thus, simultaneous observation of both DREC photons might 
contribute to the count rate in the RDEC range. The probability that both DREC photons will 
be registered as one event is estimated to be comparable with the probability of registration 
of a true RDEC event (Table l5\4j) . As the angular correlations of these two emitted photons 
have so far neither been measured nor calculated, it was assumed that the angular distribution 
of the two DREC photons is the same as for the REC photons (~ sin 2 9) and the emission 
angles for the photons are independent of each other. This assumption may cause a significant 
overestimation of the corresponding probability given in Table 15.41 It is more likely that the 
photons are emitted in opposite directions (momentum conservation) and the corresponding 
count rate is much lower than estimated above. 

For an REC cross section of 512 b the probability that one ion undergoes two collisions 
in the target, with the REC process occurring in each of them, and that both photons are 
registered as one event, is about 5.2 • 10 -14 . This, together with the mean beam current value, 
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results in a count rate of double REC photons as a single event on the order of 10 9 events/s, 
which is negligible when compared to the expected count rate of the RDEC photons (Ta- 
ble ED- 

Another possibility of simultaneous detection of two REC photons is the observation from 
two incoming ions successively undergoing REC processes in the target. This random event 
may be registered only in the time window of the order of the charge collection time in the 
silicon detector. For a typical Si(Li) crystal the collection time does not exceed 100 ns |Lip 00 . 
In case of an ORTEC silicon detector cooled with LN2 it is on the order of 10 ns |ORTb] . 
Therefore, the corresponding count rate for the mean beam intensity value is on the order of 
10~ 10 events/s and can be neglected. This estimation can be tested by the analysis of the 
Nrdec/Nrec ratio as a function of beam intensity I. The count rate of this random event 
scales with I 2 , while all the other count rates (Table lo~4"]) . including that for single REC, scale 
with I. Thus, if only true (RDEC) events are registered, the number of counts in both REC 
and RDEC regions will be proportional to the beam intensity and the ratio N re EC /^ REC 
should be independent of /. If a pile-up effect, due to sequential REC collisions was observed 
in the RDEC region, Nreec would be proportional to the square of the beam intensity and 
the N rdec /N rec ratio would increase with the beam intensity. In order to separate the data 
collected with different beam intensities registered during data taking, data files were sorted 
according to the average beam intensity and separated into three groups: 

• data taken with average beam current not exceeding 2.5 • 10 4 ions/s (average of all files: 
how = 2.0 • 10 4 ions/s) - further referred to as low intensity beam; 

• data taken with average beam current greater than 4.0 • 10 4 ions/s (average of all files: 
Ihigh = 4.3 • 10 ions/s) - further referred to as high intensity beam; 

• all the remaining data (average beam intensity: I m edium = 3.0 • 10 4 ions/s). 

This gives a beam intensity ratio Ihigh/ how = 2.15, which, if the pile-up effect was observed 
in the RDEC region, would increase the RDEC to REC ratio by a factor of almost 5. 

For the discussed beam intensities the Nrdec/Nrec ratio is shown in Fig. l5.5l for the q—1 
coincidence channel. It can be seen in this figure that the ratio does not change within the 
uncertainty limits, which means that no pile-up events were registered in the RDEC region. 
These results do not exclude the contribution of the DREC process. However, the upper limit 
of this contribution is known and will be used in Chapter [7] for further discussion. 
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Table 5.5: Calculated positions of the RDEC and REC peaks in the x-ray spectrum corresponding 
to different combinations of the initial and final states of the captured electrons. All values are given 
in keV. 







Captured target electrons 




Process — 


Valence 


One K-shell and one valence 


K-shell 


REC 


2.16 




1.88 


RDEC Is 2 


4.18 


3.91 


3.64 


RDEC lsW 


3.58 


3.31 


3.04 



5.5 Single spectra analysis 

In order to eliminate the Al K-a line from the 8+ + C data, single spectra obtained for 
bare ions without the carbon target were normalized and subtracted from the data taken with 
the foil. As the shape of the Al K-a line, despite a good statistics, was not smooth enough, the 
line in the spectrum taken without the target foil could not be normalized to the one in the 
8+ + C spectrum by simple comparison of the maxima of the lines. To establish the peak 
height of the Al K-a line in both spectra, the numbers of counts around the centroid of the 
Al K-a line in each x-ray spectrum were integrated and divided by the number of channels in 
the region of integration. This procedure prevented the contribution of the Al K-a line from 
being overestimated. Al K-a line normalized to the spectrum taken with the foil is shown in 
Fig. EH 

After subtraction of the Al K-a line, the REC structure in the single spectrum was clearly 
resolved, as shown in Fig. 15.71 For the investigated low energy collisions the binding energy 
of the target electron cannot be neglected, as it can significantly contribute to the REC photon 
energy (see Eq. 12 . 2 1) . For the considered collision system the captured electron could come 
either from the target K-shell or from the target valence band. Thus, a double peak structure 
in the REC region was expected, with the peaks separated in the spectrum by the difference 
in the binding energy of the K-shell and valence electrons in the carbon foil, which is about 
280 eV |Tho Olj . Expected positions of the corresponding REC peaks are given in Table 15.51 
As can be seen in Fig. 15.71 the obtained x-ray spectrum is well described by the sum of the 
target K-shell and valence band Compton profiles. 



44 



Anna Simon 



Correlated radiative electron capture in ion-atom collisions 



Table 5.6: Results of a x 2 test of the RDEC range of the coincidence spectra. 





N 


b 


T 


1 — OLT 


Nrdec 


P 


9-1 


567(23) 


241(13) 


25000 


<0.0001 


326(27) 


<0.0001 


q-2 


31(6) 


0.035(3) 


420 


<0.0001 


31(6) 


<0.0001 



The total number Nrec of the REC photons collected during the experiment was Nrec = 
39800(200). The uncertainty includes both the statistical error and the uncertainty generated 
during the subtraction of the background contribution. 

5.6 Coincidence spectra analysis 

In order to establish the significance of the observed structure within the RDEC energy 
range of the coincidence spectra (Fig. 15. 2} the method described in Appendix lAl was applied. 
The first step was to check whether the observed structure was due to statistical fluctuations 
(Hq hypothesis). Similar to the PIXE analysis, a = 0.05 was assumed. Then, the background 
contribution in each spectrum was obtained by integration of the REC Compton profile fit 
over the RDEC energy range (2.8-4.2 keV). This gave the total number of the background 
counts in the q - 1 and q - 2 spectra of b^ 1 = 241(13) and = 0.035(3). The total 

number of counts accumulated in the RDEC energy range during the experiment was N q ~ l = 
567(23) and N q ~ 2 = 31(6) for the single and double charge exchange channels, respectively. 
These gave the value of the statistical variable (Eq. lA~4l) of T^ 1 « 25000 and T q ~ 2 « 420, 
which corresponds to the value of olt > 0.9999 for both q — 1 and q — 2 channels. Thus, 
in both cases the hypothesis Ho has to be rejected, as the probability that the observed 
structure is due to statistical fluctuations is less than 0.0001. This means that the observed 
structure is a result of a physical process (Hi hypothesis) and the number of counts associated 
with this process is: Nj^ EC = N q ~ % — b q ~ l , where i = 1, 2, with the uncertainty given 
by AN^ EC = y^AA^-*) 2 + (A6<?-*) 2 . This leads to numbers of the RDEC counts of 
Nrdec = 326(27) and N%£ EC = 31(6) for the single and double charge exchange channels, 
respectively. Additionally, as given in Eq. IA.16| the type two error (/3) can be estimated. 
In the case discussed here, the value of /3 is less than 0.0001 which is beyond the statistical 
significance. Values of all the parameters for both single and double charge exchange channels 
are given in Table 15.61 



45 



Anna Simon 



Correlated radiative electron capture in ion-atom collisions 



Subtraction of the Compton profile based on Ref. |Big 75| from the coincidence spectra 
(Fig. 15. 2p revealed a complex structure of the RDEC line. The resulting spectra are presented 
in Figs 15.81 (c) and (d). The observed structure comprises at least two maxima which can 
be assigned to the RDEC process. It is not only a result of capture to the ground (Is 2 ) and 
excited (ls 1 2s 1 ) states of the projectile (see Fig. 15.81 (b)), but can also be attributed to the 
capture of either K-shell or valence target electrons. Combinations of the initial and final 
electron states accessible for the process and the resulting RDEC peaks positions are given 
in Table 15.51 The possible transitions are shown in Fig. 15.81 (a) with the positions of the 
lines representing their contributions to the RDEC x-ray spectrum (Fig. 15.81 (b)). For this 
presentation, probabilities of all the possible RDEC transitions were assumed to be equal. 
This spectrum should be compared with the data obtained experimentally in the single and 
double charge exchange channels after background subtraction, shown in Figs 15.81 (c) and (d). 
The latter fgure suggests that the transitions from the target valence band to the projectile 
K-shell are negligible. 

In order to make at least an estimate of the contributions of all the possible transitions, the 
sum of lines representing all the possible processes was fitted to the sum of the x-ray spectra 
registered in the q — 1 and q — 2 coincidence channels. The result is presented in Fig. 15.91 and 
the fitting parameters are given in Table 15.71 It can be noticed that in case of transitions 
which included at least one target valence electron the corresponding areas below the fitted 
curve are close to zero. Based on this observation, it can be deduced that the transitions from 
the valence band to the ls 1 2s 1 state, which overlap with the K-shell — > Is 2 transition, are 
also negligible. 

Thus, transitions from the valence band were practically not present in the RDEC spec- 
trum. Moreover, the ratio, Ra, of the areas under the curves representing transitions from 
the target K-shell to the excited and ground projectile states, Ra = A± s i 2s 2/Ai s 2, is equal to 
0.500(68) This value reflects the ratio of the RDEC cross sections for capture to the excited 
and ground states &R^Ec/ a RDEC- ^ ne theoretical value of this ratio can be estimated from 
Fig. 13.51 for the adiabacity parameter of £o = 0.82, which corresponds to the experimental 
conditions. The theoretical value of this cross section ratio is equal 0.7, which is close to the 
experimental result. 
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Table 5.7: Areas (A) of the shapes of the RDEC contributions fitted to the sum of q — 1 and q — 2 
spectra. The FWHM of all lines was set to 0.3 keV which is the width of the carbon Compton profile. 



Transition 


Peak position [keV] 


Area (A) 


2 x K-shell ls^s 1 


3.04 


118(11) 


1 x K-shell + 1 x valence -> ls 1 2s 1 


3.31 


0.0(1.9) 


2 x valence -> ls^s 1 or 2 x K-shell -> Is 2 


3.62 


237(11) 


1 x K-shell + 1 x valence -> Is 2 


3.91 


0.0(2.5) 


2 x valence -> Is 2 


4.18 


2.5(9.9) 
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E [keV] 

Figure 5.1: Experimental single x-ray spectra. In all spectra: solid line - 38 MeV 8+ + C. (a) dashed 
line - 8+ data taken without the carbon foil, (b) 8+ data after subtraction of the Al K-a line, (c) 
dotted line - 38 MeV 7+ + C, (d) dot-dashed line - 2.375 MeV protons on carbon. 
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Figure 5.2: X rays registered for 8+ + C collisions in coincidence with ions which captured (a) two 
electrons and (b) one electron. Solid line is the sum of the REC Compton profile and the Gaussian 
shape of the oxygen K-a line fitted to the spectrum. 
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Figure 5.3: Proton induced x-ray spectrum. Solid line: the RDEC range; dashed line: region consid- 
ered during background estimation. 
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Figure 5.4: Background structure in the single x-ray spectrum. The bremsstrahluug contribution 
includes all the relevant processes (SEB + AB + NB) discussed in the text. The spectrum is completely 
dominated by the REC structure. 
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Figure 5.5: Nrdec /Nrec ratio in the q — 1 coincidence spectra as a function of beam intensity. 
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Figure 5.6: 8+ spectrum taken without the carbon foil (red line) normalized to the data taken with 
the foil (black line). 
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Figure 5.7: Double structure of the REC line resolved after subtraction of the Al K-a line. 
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Figure 5.8: Possible RDEC transitions (a) and the structure of the produced x-ray spectrum (b) when 
equal cross sections for all the partial processes are assumed. Black line - the sum of all contributions. 
Additionally, corresponding RDEC spectra obtained experimentally in single (c) and double (d) charge 
exchange channels are presented. 
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Figure 5.9: The sum of spectra registered in single and double charge exchange channels with a fit of 
all possible combinations of the RDEC transitions. Fitting parameters are given in Table [5~71 



54 



Chapter 6 



The RDEC cross section 



6.1 Experimental value of the RDEC cross section 

As it was discussed in Section 15.61 RDEC photons appeared both in single and double 
charge exchange channels. The numbers of counts in the coincidence spectra, N RD } EC = 
326(27) and N q RD f EC = 31(6), in single and double charge exchange channel, respectively, 
gave the total number of RDEC counts Nrdec = 357(28). The total number of the REC 
counts collected during the experiment in the single x-ray spectrum was Nrec = 39800(200). 
The ratio Nrdec /Nrec of the total RDEC to REC numbers of counts obtained during the 
experiment was equal to 0.0092(6). One can assume that: 



do' RDEC 

Nrdec dQ 



=90° 



Nrec da rec 



(6.1) 



9=90° 



which, for the angular differential cross section for the REC process obtained from Eq. 12.61 
gives {do rdec / d£l)\e=90° = 0.71(5) b/sr. When the angular distribution of the RDEC pho- 
tons is assumed to be the same as for the REC photons (~ sin 2 6), one obtains the total 
RDEC cross section value ordec = 5.9(4) b. 

The registered REC photons can be emitted not only during the REC process but can 
also originate from the DREC process. Thus, the probability Prec of observation of a single 
REC photon is given by: 

.do DREC 



Prec fdffREa 



V dfl 



+ 2- 



Atld, (6.2) 

9=90° ' 



=90° d£l 

where d is the target thickness. In case of the DREC process the cross section has to be 
multiplied by two, as either of the two emitted photons can be registered by the x-ray detector. 
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Similarly, there is a certain probability that both REC photons emitted during theDREC 
process are registered as a single event (see discussion of this problem in Section f5.4p . Then, 
the probability Prd EC of registration of a photon in the RDEC energy range can be expressed 
as: 

Prdec=(^^-. ..An") Afld, (6.3) 



V dn 



d 2 a-£)jiEC 



6»i ,6*2=90° 



=90° dQ,idQ,2 

where the indices 1 and 2 indicate each of the emitted DREC photons. In this case: 

Nrdec Prdec 



(6.4) 



Nrec Prec 

and the RDEC differential cross section is estimated to be (damjEC /d^l)\e=9o° = 0.66(39) b/sr. 
When the angular distribution of the RDEC photons is again assumed to be the same as for 
REC photons (~ sin 2 #), the total RDEC cross section is equal to ctrdec — 5.5(3.2) b. The 
uncertainty comes mainly from the uncertainty of the target thickness given by the maker of 
the foil |ACF 09] . Both the obtained ord ec values are comparable and both are significantly 
greater than the total cross section based on the theoretical expectations {(Jrdec = 0.15 b 
[Net 09] ). 

When, according to the observations described in Section 15. 6| the contribution of the 
transitions from the target K-shell to the projectile excited state are taken into account, the 
cross section for RDEC to the ground state, tr^ EC , and for the capture to the excited state, 
a R S DECi can ^ e calculated from: 

( _ls 1 2s 1 

RDEC _ n 7 
ls 2 —u.l, 

a RDEC (6.5) 
^RDEC + a RDEC = a RDEC, 

which finally leads to the values of the RDEC cross sections <J l R s DEC and ctr^ec °^ 3.2(1.9) b 
and 2.3(1.3) b, for the capture to the ground and excited states, respectively. The cr R s ^ EC 
value is still a factor of 25 greater than the theoretical one. This leads to the ratio, R eX p, of 
the RDEC to REC cross sections of R eX p = (J rdec/ iJ R e c = 7.4(3.7) • 10~ 3 which is much 
greater than the theoretical value based on Nefiodov calculations |Nef 09] and the Stobbe 
formula for the REC cross section, R^ef = 2.9 • 10~ 4 . 
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6.2 Estimation of the (JrdecI^rec ratio in the nonrelativistic 
approach 

Simple calculations based on the principle of detailed balance can be used to estimate the 
cross sections ratio R nre l = a RDEc/ aREC ' 13-51 expresses this ratio by means of single 
and double photoionization cross sections. For a given projectile atomic number Z, the ratio 
(JDPijapi is almost constant and, for nonrelativistic photon energies, can be expressed as 
[Amu 75] : 



(6.6) 



a DPI _ 0.0932 
api Z 2 

As the photoionization cross section scales with the photon energy [Fan 59] : 

a PI ~ (-L) , (6.7) 

the ratio apj(2huj)/api(huj) = 1/32 and the cross sections ratio can be written as: 

Rnrel = = 0.00291F^. (6.8) 

&REC Z A 

In case of the discussed experiment Zt = 6, Z = 8 and the above formula gives, for F = 1, 
Rnrel = 2.3 ■ 10 -4 , which is close to the value obtained by Nefiodov. This suggests that the 
RNef value calculated by Nefiodov is not very far from the nonrelativistic estimation of the 
RDEC cross section based on the cross section for double photoionization and the principle 
of detailed balance. 



6.3 RDEC cross section based on the Yakhontov approach 

In two papers |Yak 96| lYak 97] Yakhontov calculated the value of the & r s dec /(Trec ratio 
for the 11.4 MeV/u Ar 18+ + C collisions. For a given value of the adiabacity parameter £, 
the value of this ratio scales with the projectile atomic number as Z~ 5 . In case of the argon 
and oxygen experiments, the value of the adiabacity parameter was comparable (£^ r = 0.85 
and £o = 0.82). Thus, it was possible to use this scaling for estimation of the value of Ry a kh 
for the oxygen experiment from the value calculated by Yakhontov for the argon experiment 
(Ryakh)- Since for these experiments: 



Ryakh _ Z, b 

T>Ar ry— 5 

^Yakh ^Ar 



= 57.66, (6.9) 
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Table 6.1: Comparison of the experimental values of the RDEC cross section and the R = 
a RDEC I l(T REC ratio with the ones obtained from various theoretical approaches . 





Experiment 


Nefiodov 
calculations 


Non-relativistic 
approach 


Yakhontov 
calculations 


Vrdec [b] 

R = a RDEC/ a REC 


3.2(1.9) 
7.4(3. 7)-l(T 3 


0.15 
2.9-10" 4 


0.14 
2.3-10- 4 


0.14 
2.M0~ 4 



one obtains the value of the RDEC to REC cross section ratio for oxygen of Ry a kh = 2.1 • 10 
which is consistent with both nonrelativistic estimations, i.e. Nefiodov calculations and the 
estimation based on the principle of detailed balance. 

The estimations of the RDEC cross section and the ^^j)Ecl aREC ra ti° presented in this 
chapter suggest that the experimentally obtained value of the RDEC cross section is much 
greater than theoretical estimations. Consequently, the value of the (T^^) EC cross section 
is much greater than the theoretical ones, estimated from the R = o'^jecI /(T REC ratio, as 
summarised in Table IfTTI However, according to Nefiodov [Nef 09j, the system chosen for the 
experiment does not fully meet the theoretical assumptions. Since Z ~ Zt, captured electrons 
cannot be treated as quasifree and theoretical calculations given in [Mi k 04a] IMik 04b] INef 05] 
might provide an underestimated cross section value. As the same assumptions for the cap- 
tured electrons are applied for the other theoretical approaches to the RDEC cross section 
that are given in this chapter, these calculations may also give underestimated values of the 
RDEC cross section. 
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Monte Carlo simulations of the x-ray spectra 

It was shown in Section 16.11 that the experimental value of the RDEC cross section is 
significantly greater than the one predicted by Nefiodov |Nef 09] . Thus, a Monte Carlo sim- 
ulation was performed in order to check the structure of the x-ray spectrum that could be 
obtained, assuming the value of the RDEC cross section predicted by Nefiodov was correct. 

A C++ code written by Swiat |Swi 99] was used. This is a simulation of x-ray spectra 
observed during ion- atom collisions. The x-ray spectra that can be simulated are due to 
radiative processes such as REC, DREC and RDEC as well as characteristic x-ray radiation. 
This code was thoroughly tested by its author for many collision systems |Swi 00] and was 
accepted as a reliable tool for simulations of x-ray spectra resulting from atomic collision 
processes. The simulation can reproduce almost any experimental conditions regarding beam, 
target and detector properties. It allows for implementation of any geometry of the x-ray 
detector. The parameters that have to be given in the input file are: beam diameter, target 
thickness, detector distance from the target center, observation angle with respect to the 
beam direction, detector shape and dimensions. Detection efficiency is provided as a list of 
data points which are interpolated for the given photon energy. A detailed description of 
the input parameters can be found in [Swi 00J. Here, the detector parameters were provided 
according to information presented in Chapter [4] and the detection efficiency was based on 
data given by ORTEC [ORTaJ (see Fig. 14. 5p . The code was originally designed to reproduce 
the conditions of gas-jet target experiments. Thus, by default, a Gaussian distribution of 
gas particles was used to describe the target density. Here, the code was adopted to the 
experimental conditions discussed in this thesis, a thin solid target positioned at 45° to the 
beam direction. The coordinate system used for simulation is shown in Fig. 17.11 for x = 0. 
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x-ray 
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Figure 7.1: Geometry of the experimental setup implemented in the Monte Carlo simulation, b w - the 
beam diameter, d t - target thickness in mm. The x-axis is perpendicular to the picture plane. 



It was assumed that the beam propagates in the z direction, while the center of the x-ray 
detector is positioned at x = and at positive values of y (compare Figs. 14.61 and 17. ip . The 
primary method for determination of the collision coordinates (x,y, z) within the target (see 
Chapter 3 of [Swi 00J) was replaced by: 



x = b w (Rand — 0.5), 
y = b w (Rand — 0.5), 



(7.1) 



dt(Rand — 0.5) + y tan ( — 



7T 



where b w is the beam diameter and dt is the target thickness, both given in mm. Rand E [0; 1] 
is provided by a random number generator. The beam transverse cross section was assumed 
to be a circle, thus the (x, y) coordinates were limited by: 

2 



2 i 2 ^ ; 

x +y < [ — 



(7.2) 
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Figure 7.2: Monte Carlo simulated x-ray spectra: (a) no RDEC included, (b) the RDEC cross section 
as given by Nefiodov, (c) (J^ DEC = 3 b and cr]f D 2 £ c = 2.1 b - cross sections values for which MC 
simulation gives the results closest to the experimental data, (d) Experimentally obtained singles 
spectrum. 

The REC and DREC processes were implemented with the cross sections of 512 b and 
48 b, respectively, as given in Table 15.31 As shown in the previous chapter, the RDEC 
structure observed in x-ray spectra consisted mainly of two lines. In order to simulate this 
double structure of the RDEC line (i.e. capture to the ground (Is 2 ) and excited (ls 1 2s 1 ) 
state), the input data included two kinds of the RDEC processes: capture to the ground 
state with the cross section of cr^j E( j and capture to the excited state with cross section of 
^RDEC = ^■^ (J rdec estimated according to [Nef 05] . As the photon angular distributions 
for neither RDEC nor DREC is known, they were assumed to be the same as for the REC 
photons, that is ~ sin 2 9. Additionally, the DREC photons were considered to be emitted 
independently of each other, as discussed in Section [5.41 

During the first step of the simulation, RDEC was not included in the input data. The 
obtained spectrum is shown in Fig. 17.21 (a). It can be noticed that only a few counts showed 



61 



Anna Simon 



Correlated radiative electron capture in ion-atom collisions 





50 




40 : 


c 




■3 





u 


30- 










L_ 




0) 


20- 


£ 


- 




- 


c 






10 : 




: 








50 1 


4— ' 


40 : 


c 




3 




O 


30- 


4— 




O 




!_ 




(U 

-Q 


20 : 


E 




3 




c 






10 : 




0- 



2.5 



(a) 



50 1 



•a 40- 

-1— < 

c 

§30 



.3 20- 

E 

c io- 



I 1 I 11 1 1 



3 3.5 4 
E [keV] 



4.5 



(C) 




^HjfliMi i ..J.Lln 



3 4 
E [keV] 



(b) 




3 3.5 4 

E [keV] 



Figure 7.3: The RDEC range of the x-ray spectra. Results of simulations: (a) no RDEC included, (b) 
the RDEC cross section as given by Nefiodov, (c) cr]^ DEC ~ 3 b and ct]ijjec = 2-1 b - cross sections 
values for which MC simulation gives the results closest to the experimental data, (d) Experimentally 
obtained single spectrum 



up in the RDEC energy range. These counts are due to simultaneous detection of both the 
DREC photons. Comparison with the results in Fig. 17.21 (d) shows clearly that DREC itself 
cannot explain the experimental results. Moreover, the result of the simulation shows that 
the contribution of a pile-up effect of the DREC photons is not significant as mentioned in 
Section 15.41 

The second step of the simulations was implementation of the RDEC process with the 



cross sections a 



Is* 

RDEC 



0.15 b |Nef 09] and a)£%L = Q.7a\ 



s 2 

'RDEC 



0.105 b |Nef 05j f or the 



capture to the ground and the excited state, respectively. Again, the results showed that this 
was not sufficient to explain the experimental data (compare Figs 17.21 (b) and (d)). 

During the third step of the simulation the RDEC cross section was increased until the 
resulting spectrum was comparable with the experimental one (Fig. 17.21 (d)). The best results, 
shown in Fig. 17.21 (c), were obtained for (y x ^ DEC = 3 b and (t r s ^^q = 0-7<7^> EC < = 2.1 b. 
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Table 7.1: Ratios of the numbers of counts in the RDEC and REC regions obtained during Monte 
Carlo simulations compared with the experimental value. 





RDEC cross section [b] 


Number of counts 


RDEC/REC ratio 




Is 2 , rr ls 1 2s 1 
"RDEC + "RDEC 


REC range 


RDEC range 




1 


40100(180) 


71(8) 


0.0018(2) 


simulation < 


0.15 + 0.105 


39900(180) 


225(8) 


0.0056(5) 




3 + 2.1 


39700(180) 


449(20) 


0.0113(3) 


experiment 


3.2 + 2.3 


39800(200) 


357(28) 


0.0092(6) 



However, it can be noticed that the RDEC structure within the experimentally obtained 
spectrum (Fig. 17.21 (d)) is placed on the REC tail, which seems to be much broader than 
the simulated one. This may be due to the fact that for the simulation the Compton profile 
of diamond was used as given by Reed |Ree 72] . which might differ from the one for the 
amorphous carbon foil that was used during the experiment. However, the main purpose of 
the simulation was to investigate the ratio of the numbers of counts related to the RDEC and 
REC processes, thus the Compton profile here is of a minor importance. 

In Fig. 17.31 only the RDEC range of all the spectra was shown on the linear scale, for 
better visualization of the double structure of the RDEC line. As a test of the result of the 
simulations, the ratio of the numbers of counts in the RDEC and REC range was calculated. 
The obtained values are shown in Table 17.11 It can be seen that the results of the third step 
of the simulation are in the best agreement with the experimental data. This again shows 
that the cross section value calculated by Nefiodov, even when capture to the excited state is 
included, is insufficient to explain the experimental results. 
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Chapter 8 



Conclusions 



In this dissertation an experiment dedicated to the radiative double electron capture pro- 
cess (RDEC) was presented. The experiment was carried out at Western Michigan Univer- 
sity using the 6 MV tandem Van de Graaff accelerator. The choice of the collision system, 
8+ + C at 38 MeV, was based on the recent theoretical calculations of the RDEC cross 
section |Mik 04a] IMik 04b"l INef 05] . which pointed to mid-Z ions and low collision energy as 
the best systems for observation of RDEC. The theoretical approach also suggested an en- 
hancement of the RDEC cross section in such systems due to capture of electrons from the 
target valence band. Moreover, it was pointed out that the capture to the excited state of the 
projectile may be a significant contribution to the process. 

The results allowed for the first experimental verification of RDEC and provided a test 
of the main theoretical predictions |Mik 04a] IMik 04b] INef 05| . The obtained x-ray spectra 
revealed a complex structure of the RDEC line. However, the capture from the target valence 
band, which, according to the theory, was supposed to significantly contribute to the RDEC 
process, was not confirmed. The observed structure allowed for identification of capture to 
the projectile ground (Is 2 ) and excited (ls 1 2s 2 ) states. The ratio of the counts which could 
be associated with these two processes gave a ratio of the RDEC cross sections for the capture 
to the excited and the ground states &r^ec / a RDEC = 0.500(68), which is close, considering 
the data statistics, to the theoretical value of 0.7. 

The latter value, together with the observed ratio of the numbers of counts in the RDEC 
and REC range of the x-ray spectra, NrdeC/ 'NreC = 0.0092(6), allowed for estimation of 
the RDEC cross sections: 

• &REIEC = 3-2(1.9) b for the capture to the ground state, 
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Table 8.1: Summary of the results of the theoretical calculations and the experiments dedicated to 
the RDEC process. 



z 


E [MeV/u] 




Zi 




Vrdec l mb l 




Ref. |Mik 04a| 


Ref. [Yak 97| 


experiment 


18 


11.4 


0.840 


6 


3.2 


1.85 


<5.2 [War 95) 


92 


297 


0.841 


18 


2.5-10~ 2 


5000 


<10 (Bed 03] 


8 


2.375 


0.820 


6 


1.5-10 2 (a) 


1.4-10 2 W 


3.2(1.9)-10 3 (c) 



( Q ) provided by Nefiodov |Nef 09] 

( b ) estimated from the Ry a kh = a R S DEC 1 1(7 REC ratio 

( c ) this experiment 



• cr)j£)^(7 = 2.3(1.3) b for the capture to the excited state. 

The main results of the experiment have already been published |Sim 10] . The obtained value 
of the cross section is a factor of about 25 greater than the theoretical value provided by Ne- 
fiodov |Nef 09] . Consequently, the ratio R = cr 1 ^ )EC /aREC = 0.0074(37) is also significantly 
greater than estimated from the theory, as shown in Table 18.11 

The results of the so far conducted experiments dedicated to the RDEC process are given in 
Table [87T1 . The discrepancies between the theories and experiments, as well as the differences 
between various theoretical approaches, show that further investigation of the RDEC process 
is necessary. A similar experiment chould be performed, with the same experimental setup 
available at WMU, but with an additional absorbing material in front of the x-ray detector. 
This will exclude all the pile-up effects and contributions of the DREC process, which formed 
a crucial problem of the background analysis in the presented experiment. Moreover, the 
angular distribution and correlations of the DREC photons should be measured as DREC 
interferes with other radiative capture processes. 

As the role of the capture from the valence band of the solid target was not verified, 
application of a gas jet target with a light gas (He) would be useful, due to a significant 
reduction of background and a much smaller probability of multiple collisions. Additionally, 
a system with a greater difference between projectile and target atomic number (for example 
Ca 20+ , Ar 18+ or Xe 54+ on He) would allow for better separation of photons originating from 
the capture to the ground and excited state of the projectile. All these conditions can be 
satisfied during an experiment on the ESR gas jet target at GSI. 



65 



Appendix A 



Statistical analysis of the observed signal 



In order to check if the observed structure within a given region of the spectrum is a 
result of a physical process or just a statistical fluctuation, a method suggested in |Ead 89] 
was applied. Here, a brief description of the approach is presented. 

Alternatively, one can calculate the probability of observation of a statistical fluctuation 
within the range of interest. This method assumes that the background distribution is known. 

A null hypothesis can be defined as: 

Hq: there is no physical effect within the AB range. 

It is assumed that the background shape can be described by a function 6(x, 0), which 
depends on the observed variable x and unknown parameters 0. In this case the number of 
the background counts in the AB region (Fig. IA.1|) is given by: 

b AB = [ 6(x,0)dx. (A.l) 

J A 

As oab is a function of estimators 0, the variance can be obtained by a substitution of 
variables: 

a\ B = D T VT>, (A.2) 

where D denotes a vector of derivatives: 

db AB <' B 



L ^ M)ix (A - 3) 



If Nab denotes the number of events in the AB range, the optimal test statistic, which checks 
if Nab significantly differs from 6ab ; is: 

m (Nab - b A B? 



V(N AB - b AB ) 



(A.4) 
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Figure A.l: Example of an experimentally obtained spectrum with a structure within the AB range. 



where V(Nab — i>A B ) is a variance. If the Hq hypothesis is true, then: 

E(N AB ) = b AB 



(A.5) 



V{N AB ) = b AB (A. 6) 

where E(N AB ) is the estimator of the number of counts in the AB region and the estimator 



of b AB is b AB . Thus: 



V{N AB - b AB ) « b AB + u\ B - 2cov(N AB , b AB ). 



(A.7) 



As during the estimation of O parameters the AB range was excluded, N AB and b AB are 
uncorr elated: 

V(N AB - b AB ) w b AB + a\ B (A.8) 

and 



T = 



(Nab - b AB f 
o\ B + b AB 



(A.9) 



If the number of counts N AB is large, it has a normal distribution around b AB and the statistic 
T behaves as a x 2 distribution for DoF = 1 (degrees of freedom). 



67 



Anna Simon 



Correlated radiative electron capture in ion-atom collisions 



Table A.l: Quantiles of the x 2 distribution for DoF = 1 [Karnj . 



a 0.005 0.010 0.020 0.025 0.050 0.100 0.250 0.300 0.500 

X 2 4.0E-5 0.0002 0.0006 0.0010 0.0039 0.0642 0.1015 0.1485 0.4549 

a 0.700 0.750 0.800 0.900 0.950 0.975 0.990 0.995 0.999 

y 2 1.0742 1.3233 1.6424 2.7055 3.8415 5.4119 6.6349 7.8794 10.8270 



For the obtained value of the statistical variable T the value of olt can be estimated 
based on Table IA.1I If (1 — olt) > ol, where a is the probability of the type one error, the 
hypothesis can be accepted. The deviation is small enough for chance alone to account for it. 
If (1 — olt) < ol, the hypothesis should be rejected, as there may be some factor other than 
chance operating for the deviation to be so great. 

If, based on statistics IA.9| Hq was rejected, estimation of the signal s = Nab — bAB can 
still be used, but with a different estimator variance. The problem one is about to solve is 
now a test of hypothesis H\ against hypothesis Hq. 

Hi: there is a physical signal s and background bAB in the AB range. 

In this case the variance will be given by: 



where, for the same reason as in Eq. IA.8| covariance is not included. The risk of making a 
type two error is equal to: 



V(N AB - b AB ) « N AB + o\ B , 



(A.10) 



P = P(d<X a \H 1 ), 



(A.ll) 



where d = VT and X a is defined by 3>(A Q ) 



a, $ being a cumulative normal distribution 



function, given by: 




(A.12) 



For the H\ hypothesis the mean value and variance are given by: 



E(Nab) = V{N AB ) = b AB + s 



(A.13) 



while d has a normal distribution N(/x, a 2 ), where: 



s 



(A.14) 
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a 



bAB + s + a 



bAB + v\ B 



AB 



Thus, finally: 



= $ 



A«\/ bAB + a 2 



AB 



bAB + S + <T^ B 



(A.15) 



(A.16) 
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